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ABSTRACT 

The NGC 1333 region was observed in the H2 1-0 S(l) line. The H2 images cover a 5' x 7' region around IRAS 4. 
Numerous H2 emission features were detected. The northeast-southwest bipolar outflow driven by IRAS 4A was studied 
by combining the H2 images with SiO maps published previously. The SiO-H2 outflows are continuous on the south- 
western side but show a gap on the northeastern side. The southwestern outflow lobe curves smoothly, and the position 
angle increases with the distance from the driving source. The base and the outer tip of the northeastern outflow lobe 
are located at positions opposite to the corresponding parts of the southwestern lobe. This point-symmetry suggests that 
the outflow axis may be drifting or precessing clockwise in the plane of the sky and that the cause of the axis drift may 
be intrinsic to the outflow engine. The axis drift model is supported by the asymmetric lateral intensity profile of the 
SiO outflow. The axis drift rate is ^0?01 1 yr^^. The middle part of the northeastern outflow does not exactly follow the 
point symmetry because of the superposition of two different kinds of directional variability: the axis drift of the driving 
source and the deflection by a dense core. The axis drift model provides a good explanation for the large deflection angle 
of the northeastern outflow. Other H2 emission features around the IRAS 4 region are discussed briefly. Some of them 
are newly found outflows, and some are associated with outflows already known before. 

Subject headings: ISM: individual (NGC 1333 IRAS 4) — ISM: jets and outflows — ISM: structure — stars: formation 



1. INTRODUCTION 

The NGC 1333 region contains numerous young stellar ob- 
jects and outflows (Aspin et al. 1994; Hodapp & Ladd 1995; 
Rodriguez et al. 1999). It is one of the most active sites of star 
formation in the solar neighborhood, and the large number of 
active outflow sources was described as a "microburst" of star 
formation (Bally et al. 1996). Several protostars are located in 
the NGC 1333 molecular cloud, and IRAS 4A is one of the 
brightest submillimeter sources among them (Sandell & Knee 
2001). IRAS 4A is a Class protobinary system (Sandefl et al. 
1991; Lay et al. 1995; Looney et al. 2000). 

Single-dish observations showed an interesting directional 
variability of the IRAS 4A molecular outflow (Blake et al. 
1995). Detailed structure of the IRAS 4 A molecular outflows 
within ^1' of the driving source was revealed by interferomet- 
ric observations. The interferometric maps showed that the ap- 
parently large change of outflow direction is partly owing to the 
existence of two outflows driven by IRAS 4A (Choi 2001b) and 
partly owing to the sharp bend (or deflection) in the northeastern 
lobe of the main bipolar outflow (Girart et al. 1999; Choi 2005a, 
hereafter Paper I). Choi (2005a) suggested that the sharp bend 
was caused by a collision between the northeastern outflow and 
a dense core in the ambient molecular cloud. The jet-core col- 
Usion model was supported by several lines of evidence includ- 
ing the asymmetric morphology of the main outflow, the HCN 
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emission from shocked high-density gas near the impact point, 
and the high-density core located in the course of the north- 
eastern outflow (Paper I; Choi 2005b). However, the deflection 
cannot completely explain the difference of position angles be- 
tween the SiO outflow near the driving source and the large- 
scale (> l')CO/H2 outflow. 

To investigate the overall structure of the NGC 1333 IRAS 
4A outflow system, we observed the outflows in several molec- 
ular tracers. In Paper I, we presented high angular resolution 
observations in the SiO i; = 0J=1^0 line using the Very 
Large Array (VLA). In this paper, we present our observations 
in the H2 1-0 S'(l) line using the University of Hawaii (UH) 2.2 
m telescope and the Subaru telescope. In § 2 we describe our 
H2 observations. In § 3 we report the results of the H2 imaging. 
In § 4 we compare the H2 and the SiO images and discuss the 
variability of the IRAS 4A main outflow. 

2. OBSERVATIONS 

2.1. UH 2.2 m Telescope Observations 

The molecular hydrogen emission features associated with 
IRAS 4 were included in a wide field image of the NGC 1333 
region in the H2 1-0 5'(1) emission line obtained in the nights 
of 1996 July 30 to August 3 using the QUIRC camera (Hodapp 
et al. 1996) at the UH 2.2 m telescope. The individual integra- 
tion time was 300 s. The image was produced by stepping the 
telescope in steps of 30" in the north-south direction. There- 
fore, most positions in the image were covered by six individ- 
ual frames for a total integration time of 30 min. In the overlap 
region between adjacent vertical stripes, the coverage is 12-fold 
for a total integration time of 60 min. 

The UH H2 image shows a set of artifacts that are the re- 
sult of the residual excess dark current of the HgCdTe array 
used for the observations. This particular array, one of the ear- 
liest HAWAII-I devices, showed a strong excess of dark current 
in pixels previously exposed to high signal levels, with a de- 
cay time of the order of many minutes. As a result of this ef- 
fect, combined with our method of collecting the individual im- 
ages of this mosaic in a north-to-south stepping pattern, bright 



sources, most prominently SVS 13, are associated with a long 
string of residual images south of them, in 30" spacing, and 
with slowly decaying intensity. The faint extended nebulosity 
discussed in this paper is not affected by this artifact. 

2.2. Subaru Telescope Observations 

Imaging was carried out on 2000 September 13 using the 
Cooled Infrared Spectrograph and Camera for OH-Airglow 
Suppressor (CISCO), which was equipped with a 1024 x 1024 
HgCdTe detector covering a field of view of 1 10" x 1 10" with 
a plate scale of O'.'U pixel"^ (Motohara et al. 2002). We ob- 
tained narrow band images with and without H2 1-0 S(l) line 
emission. The H2 1-0 filter was centered on A = 2. 120 /im with 
a width of 0.020 fim (FWHM), while the N215 filter for contin- 
uum emission was centered on 2.147 /im with a width of 0.021 
/xm (FWHM). Sixteen exposures with 60 s each were made for 
each of the filters by stepping the telescope in steps of 10"-20" 
along the east-west and/or north-south directions, resulting the 
total on-source exposure time of 960 s for the central 90" x 
90" region. Each frame of H2 1-0 and N215 filter images was 
flat-fielded using sky-flat data and was subtracted by a corre- 
sponding sky frame, which was made for each of the two filters 
with median filtering. It was found that all the extended fea- 
tures in the H2 1-0 filter image arise from the H2 1-0 S(l) line 
by comparing it with the N215 filter image. We thus use only 
the H2 1-0 filter image in this paper. 

2.3. Astrometric Calibration 

Accurate astrometry is essential to comparing the H2 image 
with the SiO image. Since the adjacent continuum emission was 
not subtracted, the H2 images include both the line and the con- 
tinuum emission. In the UH H2 image, 20 field stars listed in the 
2MASS Point Source Catalogue** were identified. In the Subaru 
H2 image, only three 2MASS stars were found. 

With the 2MASS coordinate data, the geomap package of 
IRAF was used to find the spatial transformation functions. 
Then the geotran package was used to transform the observed 
images. Comparing the stellar coordinates of the transformed 
UH H2 image and of the 2MASS catalogue, the rms of position 
difference is 0'.'19. This value is about the same order of the 
position uncertainties in the 2MASS catalogue, which confirms 
that the transformation is acceptable. 

Angular resolutions of the resulting images were measured 
by fitting Gaussian intensity profiles to stellar objects. The an- 
gular resolutions (FWHM) are 0'.'82 for the UH H2 image and 
(y.'46 for the Subaru H2 image. 

2.4. VIA Data 

Details of the VLA observations and the results were pre- 
sented in Paper I. The NGC 1333 IRAS 4 region was observed 
using VLA in the SiO v — Q J — 1 ^r {) line. The resulting im- 
ages have a restoring beam of FWHM = 1"96. The SiO image 
covers a region up to ^47" from IRAS 4A. 

3. RESULTS 

Figure 1 shows the UH H2 image of the NGC 1333 cloud 
covering a 5.'2 x 7.'5 region. The Subaru H2 image (Fig. 2) has 

'This publication makes use of data products from the Two Micron All Sky Sur- 
vey (2MASS), which is a joint project of the University of Massachusetts and 
the Infrared Processing and Analysis Center/California Institute of Technology, 
funded by the National Aeronautics and Space Administration and the National 
Science Foundation. 



a better sensitivity and a higher angular resolution than the UH 
image, but the field of view is smaller In Figures 36 and 4a, the 
VLA SiO map of the IRAS 4A outflows is superposed on the 
H2 images for comparison. 

The structure of the H2 emission features in the NGC 1333 
region was previously discussed by Hodapp & Ladd (1995). 
Recent observations in the submillimeter and centimeter con- 
tinuum (Rodriguez et al. 1999; Sandell & Knee 2001) provide 
helpful information for understanding the relation between the 
H2 emission features and young stellar objects in this region. 
The NGC 1333 region is crowded with outflows, and the H2 
emission features should be inspected carefully to avoid confu- 
sion. 

Since the main interest of this paper is the structure of the 
IRAS 4A molecular outflows, we describe the H2 emission fea- 
tures related to these flows first. While the H2 images show the 
large-scale structure of the main outflow, little H2 emission can 
be seen near the driving source, probably owing to a high ex- 
tinction through the dense core. In contrast, the SiO line nicely 
shows the outflow structure near the driving source, probably 
owing to a high density. Therefore, the H2 images and the SiO 
map provide complementary information. 

On the northeastern side, there is a relatively wide (^35") 
gap between ASR 57 and the northeastern end of the SiO out- 
flow (Fig. 3&). This gap is probably real because CO outflow 
maps also show a similar gap (see Fig. 3 of Blake et al. 1995). A 
chain of H2 features (ASR 57, HL 10, and HL 1 1 = HH 347 A) 
seems to belong to the IRAS 4A main outflow. Each of these 
features is elongated in the general direction of the flow. HH 
347B is located away from the chain of H2 features, and it prob- 
ably does not belong to the IRAS 4A outflow. 

On the southwestern side, the blueshifted SiO outflow is con- 
nected with the H2 outflow, and there is a considerable overlap 
(Fig. 4a). Near the driving source, the faint H2 feature CHH 
15 coincides with SiO outflow peak (SiOOP) 7. The H2 feature 
CHH 17 coincides with SiOOP 10 and 11 within a few arc- 
seconds, and the northeastern end of CHH 18 coincides with 
SiOOP 12 (Fig. 4a). Further away from the driving source, a 
chain of H2 features (HL 5/3) seems to belong to the IRAS 4A 
main outflow. 

Along the southern SiO outflow, the H2 feature HL 6 is lo- 
cated in the middle of SiOOP 14, 15, and 16. Unlike the main 
outflow, the southern outflow shows no H2 emission feature be- 
yond the SiO outflow lobe. Therefore, the VLA SiO map prob- 
ably traces the whole extent of the southern outflow. 

3.1. Proper Motion 

Proper motions of H2 features were measured by compar- 
ing the Subaru and the UH images. The coordinate systems 
of the two images were aligned by comparing the positions of 
three stars near the H2 features (ASR 53, ASR 103, and CHH 
18 IRS). Since these stars are visible only in the near-IR or 
longer wavelengths, they are either stars embedded in the NGC 
1333 molecular cloud or background stars, and their proper 
motions may be negligible. For bright H2 features with well- 
defined peaks, position offsets were measured by comparing 
pixels around the peaks, typically pixels within ^2" from the 
peak positions. Additional measurements were made by com- 
paring the Subaru image with the H2 image of Hodapp & Ladd 
(1995) in a similar manner, which provides a longer time base- 
line but a larger uncertainty owing to the worse angular reso- 
lution. The two sets of measurements were consistent to each 
other, and the average proper-motion vectors are shown in Fig- 
ure 5. The uncertainty of position offset was measured from the 
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Fig. 1. — Map of the H2 1-0 5(1) line toward the NGC 1333 IRAS 4 region from the UH 2.2 m telescope. Extended sources are labeled. HL 
source numbers are from Hodapp & Ladd (1995). ASR source numbers are from Aspin et al. (1994), and Herbig-Haro (HH) object numbers are 
from Bally et al. (1996) and references therein. Newly identified sources are numbered with a prefix CHH. Squares: Centimeter continuum sources 
labeled following Table 1 of Rodriguez et al. (1999). VLA 25 is associated with IRAS 4A, and VLA 4 is associated with SVS 13. Open circles: 
Submillimeter continuum sources labeled following Table 1 of Sandell & Knee (2001). 

A preprint containing high-quality images of Figs. 1, 2, 36, and 4a can be found at http : //hanul . kasi . re . kr/"minho/ Publications .html. 
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Fig. 2. — Map of theH2 1-0 S(l) line toward the NGC 1333 IRAS 4 region from the Subaru telescope. Extended sources are labeled. Filled circles: 
Millimeter continuum sources in the IRAS 4 region (Looney et al. 2000). 



residual position differences of the five stars common to both 
images (the three stars above and ASR 54/55). 

The proper-motion vectors of HL 4 peaks show a system- 
atic motion toward southeast, and those of HL 5 peaks show a 
motion toward southwest. Assuming a distance of 320 pc (de 
Zeeuw et al. 1999), the proper motion of HL 5 peaks ranges 
from 30 to 120 km s~^ with an uncertainty of 50 km s^^. In 
the discussions below (§ 4), a proper motion of 100 km s^^ or 
O'.'O? yr^^ is assumed for calculations of the timescale of the 
IRAS 4A main outflow. 



4. DISCUSSION 

4.1. Large-Scale Variability of the IRAS 4A Main Outflow 

As discussed in Paper I, the northeastern SiO outflow shows 
a directional variability that is probably caused by a collision 
with a dense core. On top of this variability, yet another direc- 
tional variability is revealed by combining the SiO map with 
the H2 images. This variability appears to affect both the south- 
western and the northeastern outflows in a large (^5') scale. 
Since the northeastern outflow is more complicated owing to 
the superposition of the two different variabilities, we describe 
the southwestern outflow first. 
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Fig. 3. — (a) Position angle of the emission peaks in the IRAS 4A main outflow plotted as a function of the angular distance with respect to 
the driving source (IRAS 4A2). Top panel: Southwestern outflow peaks. Bottom panel: Northeastern outflow peaks. Crosses: SiO outflow peaks 
(Paper I; see Fig. 4). The vertical size corresponds to the beam (FWHM). Open circles: H2 emission peaks identified in the Subaru image. Filled 
circles: H2 emission peaks identified in the UH image. Dotted straight line: Linear least-squares fit to the southwestern outflow peaks. Dashed 
cur-ve: Parabolic least-squares fit to the southwestern outflow peaks. In the bottom panel, the fits in the top panel were drawn in reverse, i.e., with the 
position angle shifted by 1 80° . (b) Maps of the H2 line (heat scale) and the SiO line (contours) toward the NGC 1333 IRAS 4 region. The SiO map 
is the same as the one shown in Fig. 3a of Paper I. Dashed cun'e: The parabolic fit in (a) drawn in the image coordinate. Large cross: HCO+ peak 
position of the dense core obstructing the northeastern outflow (Choi 2005b). Filled circles: Millimeter continuum sources (Looney et al. 2000). 



The southwestern lobe of the SiO outflow appears very 
straight: all the SiO outflow peaks are located within a beam size 
from a straight line (see Fig. 2a of Paper I). When the SiO and 
the H2 images are superposed, however, the southwestern out- 
flow appears to curve smoothly. That is, the position angle (PA) 
of the outflow peaks with respect to the driving source (IRAS 
4A2) increases with the angular distance (Fig. 3a top panel). 
Polynomial fits were drawn in Figure 3 to outline the flow. A 
change of flow direction can be caused either by intrinsic vari- 
ability, such as geometric change of the outflow engine, or by 
external perturbation, such as density gradients in the external 
medium. The nearly constant gradient of the PA-distance rela- 
tion suggests that the large-scale variability is probably owing 
to an intrinsic variability. Another important clue to the vari- 
ability mechanism may come from the symmetry of the outflow 
morphology. 

Near the driving source, the direction of the northeastern out- 
flow is exactly opposite to that of the southwestern outflow. 
The northeastern outflow changes direction abruptly at a point 
^23" away from the driving source, and the position angle in- 
creases at a very high rate (Fig. 3 bottom panel). At a large dis- 



tance, however, the PA-distance relation converges back to the 
reverse of the southwestern flow. Therefore, it appears that the 
southwestern and the northeastern outflows would have shown 
a change of flow direction in a point-symmetric way, if it were 
not for the jet-core collision of the northeastern outflow. The 
symmetric morphology suggests that the change of flow direc- 
tion in the large scale is caused by a geometric reconfiguration 
of the outflow engine. That is, the axis of the IRAS 4A2 accre- 
tion disk may be changing direction, and the outflow axis may 
be drifting'' together The axis appears to drift clockwise in the 
plane of the sky. 

The IRAS 4A outflow is not the only example of point- 
symmetric, or S-shaped, outflow driven by young stellar ob- 
jects. Other examples include the HH 340/343 outflow driven 
by NGC 1333 J032845.3H-3 10542 (K) and the H2 jet driven by 
V380 Ori NE (Hodapp et al. 2005; Davis et al. 2000). Some 
outflows driven by high-mass (proto)stars appear to have simi- 

"Outflows with changing directions are often described as "precessing". In this 
paper, we use the term "drifting" instead of precessing because there is no evi- 
dence that the change is periodic. 
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lar shapes, e.g., IRAS 20126H-4104 (Shepherd et al. 2000). 

4.2. Lateral Structure of the SiO Outflow 

The lateral structure of the SiO outflow was briefly described 
in § 3.3 of Paper I. Here we describe the structure in more detail 
and discuss it in the context of the drifting axis. Figure Ah shows 
the lateral intensity profiles. The SiO intensity profiles of the 
undeflected part of the main outflow (peaks 4-12) usually show 
a steep slope on the clockwise side and a relatively slower slope 
(or shoulder) on the counterclockwise side. 

The asymmetric slopes of the intensity profiles can be ex- 
plained by the drift of the flow axis. Since the axis of the pri- 
mary jet drifts clockwise, the ambient molecular gas on the 
leading (clockwise) side would be relatively fresh (and proba- 
bly denser), while the molecular gas on the opposite side would 
have been already disturbed and accelerated in the past. The 
shear zone between the primary jet and the ambient medium 
may be thinner on the leading side than on the trailing side. 
Then the SiO molecules in the shear zone of the leading side 
would be more strongly shocked, excited more highly, and make 



stronger emission, which could produce the asymmetric inten- 
sity profile across the outflow lobe. On the trailing side, the SiO 
emitting gas may be gradually cooling down as the jet drifts 
away. Therefore, the asymmetric lateral intensity profile pro- 
vides additional support for the drifting axis model. 

4.3. Drifting Axis of the IRAS 4A Main Outflow 

The PA-distance relation of the main outflow (Fig. 3a) sug- 
gests that the amount of axis drift within the observed extent of 
the main outflow (i.e., from the vicinity of the driving source to 
HL 3) is --20° and that the rate of drift is --10° arcminute^^ 
Assuming a proper motion of 0'.'07 yr^^ (see §3.1), the tempo- 
ral rate of axis drift is --0?01 1 yr^^. 

Several authors suggested possible mechanisms of intrinsic 
outflow variability (e.g., Eisloffel & Mundt 1997; Fendt & Zin- 
necker 1998; Shepherd et al. 2000), but most of them cannot 
generate a large amount of drift as seen in the IRAS 4A out- 
flow. The tidal interaction between IRAS 4A1 and A2 cannot 
explain the drift because the separation between them (> 540 
AU) is too large to cause a strong effect. In principle, A2 itself 
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could be a close binary system that has not been observation- 
ally resolved yet. Assuming a proper motion of (X'O? yr~^, the 
130" separation between IRAS 4A and HL 3b gives an outflow 
timescale of '^2000 yr. Since the flow is not showing multiple 
turns, the precession period would be longer than 4000 yr, but 
it should be shorter than ^16000 yr to keep the precession an- 
gle reasonably small, say less than 90°. Since the precession 
period would be longer than the binary orbital period by a fac- 
tor of '^20 (Bate et al. 2000), the expected orbital period of the 
unresolved binary system is 200-800 yr. The corresponding bi- 
nary separation would be 30-80 AU, assuming that the mass of 
IRAS 4A2 is -0.9 Mq, which is about half the mass of IRAS 
4A (mass from Choi 2001b, scaled to 320 pc). Then the angular 
separation would be 0'.'09-0'.'25, which is too small for exist- 
ing millimeter interferometers to resolve. The accretion disk in 
such a system, however, would probably be too small to drive 
an energetic outflow. 

Another possible mechanism is the anisotropic accretion that 
can occur when the angular momentum of currently accreting 
mass is not parallel to the rotation axis of previously accreted 
matter. Considering that IRAS 4A is a binary system, we pro- 
pose the following scenario. The initial IRAS 4A core had frag- 
mented into two subcores (Al and A2), and each of them col- 
lapsed to form a protostar and started an outflow activity. While 
the total angular momentum of the A1/A2 system may be par- 
allel to the rotation axis of the initial core, the spin axis of each 
individual protostar (and the initial direction of each outflow) 
could have been oriented differently. (This scenario is consis- 
tent with the observations because there is a considerable dif- 
ference between the directions of the main outflow driven by 
A2 and the southern outflow driven by Al.) As the protostar ac- 
cumulates mass from the large outer envelope, the axis of the 
accretion disk would shift gradually and converge to the direc- 
tion of the average angular momentum of the whole protostellar 
core. This scenario can be tested in principle by observation- 
ally determining the rotation axis of the IRAS 4A core, but the 
velocity structure of this region may be too complicated (Choi 
et al. 2004) for such a test. In the anisotropic accretion model. 



the change of direction may happen only once, i.e., may not 
be periodic or oscillating, which can explain why most of the 
direction-changing outflows observed so far do not show multi- 
ple turns. 

Yet another possibility is the misalignment between the ro- 
tation and magnetic fields. Numerical simulations show that the 
rotation axis, magnetic fields, and disk orientation can precess 
significantly (Matsumoto & Tomisaka 2004). In such systems, 
changes of outflow direction may be observable. The evolution 
of the field direction with respect to the rotation axis may de- 
pend on the relative strength between the rotation and the mag- 
netic fields (Machida et al. 2006). 

4.4. Deflection of the Northeastern Outflow 

Theoretical models of jet-cloud collision predict that, in the 
simple cases of steady jet, the deflection process finishes in a 
few hundred years either because the jet penetrates the dense 
clump or because the jet is pinched off (Raga & Canto 1995; 
de Gouveia Dal Pino 1999). When the whole extent of the H2 
outflow is considered, the timescale of the IRAS 4A main out- 
flow is ^^2000 yr. Then the deflection process of a steady jet 
would have been almost finished in this timescale, and the de- 
flection angle would have been small. The large deflection an- 
gle observed in the northeastern outflow suggests that the real 
situation is more complicated than the simple models. It was 
suggested that the deflection timescale can be longer if the jet 
beam is not steady so that the impact point does not stay at one 
position (Raga et al. 2002). Therefore, the drifting axis model 
provides a good explanation for the large deflection angle cur- 
rently observed in the IRAS 4A outflow. 

Assuming that the width of the (undeflected) primary jet 
at the impact point is smaller than 4'.' 3 (deconvolved width of 
SiOOP 4), and considering the axis drift rate in the previous 
section, it takes less than —1000 yr for the impact point at the 
"surface" of the dense core to move across a jet width. Since this 
crossing time is comparable to or shorter than the cloud pene- 
tration timescale, the large deflection angle can be maintained. 

The deflection history of the northeastern outflow can be 
summarized in the following scenario. Choi (2005b) suggested 
that the dense core obstructing the northeastern outflow is lo- 
cated ~50" from IRAS 4A2 with a position angle of —16° 
(Figs. 36 and 6). The detection of HCO+ and HCN emission 
suggests that the obstructing core is very dense (—10^ cm^'^) 
(Paper I; Choi 2005b). In the distant past, the direction of the 
outflow had a large position angle, and the flow would have 
passed by the eastern side of the obstructing core with a large 
impact parameter and would have been deflected only slightly. 
For example, HL 1 1 would have had an (projected) impact pa- 
rameter of —19" (or 6000 AU at 320 pc). As the outflow axis 
drifts clockwise, the impact point would move on the "surface" 
of the dense core to the west, and the impact parameter de- 
creases. Consequently, the outflow would be deflected more and 
more severely, and the deflection angle may increase with time. 
The position angle of the outflow axis in the vicinity of the driv- 
ing source is —19°, and the jet currently colliding against the 
dense core may have an impact parameter of —3" or 900 AU. 
That is, the northeastern outflow is currently making a nearly 
head-on collision against the obstructing core, and consequently 
the deflection angle is large. 

Reipurth et al. (1996) suggested that the HH 270/1 10 system 
may be an example of deflected outflow. This interpretation, 
however, was questioned because high-density molecular gas 
was not detected near the proposed impact point (Choi 2001a). 
Raga et al. (2002) suggested that the apparent conflict can be 







RA offset (arcsec) 

Fig. 6. — Maps of the H2 line (gray scale), the SiO line (thick contours. 
Paper I), and the 8 km s^^ component of the HCO"'" J — 1 ^ line 
(thin contours, Choi 2005b). The HCO+ and the SiO maps are not 
corrected for the primary beam response. Shown at the bottom left- 
hand comer is the synthesized beam of the HCO^ map. The markers 
are the same as those in Fig. 36. 



reconciled if the impact region is located far away from the 
crossing point of the HH 270 and the HH 1 10 axes. If the HH 
110 flow is indeed a deflected segment of the HH 270 flow, its 
deflection history must be more complicated than the case of 
IRAS 4A. More detailed and sensitive observations are needed 
to understand the HH 270/1 10 system better. 



CHH 7/10/14. — These are isolated H2 features with fuzzy 
appearances. Their nature is unclear. CHH 10 is associated with 
HH 759 (Walawender et al. 2005). 

CHH 8/9. — These H2 features seem to belong to the HH 5 
system. 

CHH 11. — This is a system of weak H2 features. It appears 
to be related with HL 8. 

CHH 12/13. — These H2 features probably belong to a single 
outflow. A probable driving source is SK 1 . 

CHH 16. — This is a faint H2 feature associated with SiOOP 
17. Since SiOOP 17 is kinematicafly distinct from the IRAS 4A 
main outflow, it is not clear if it is related with the IRAS 4A 
system. See § 3.1 of Paper I for details. 

CHH 1 9.— This H2 feature probably belongs to the IRAS 4A 
outflow. 

HL &— Hodapp & Ladd (1995) suggested that HL 8 may 
be an extension of HL 4. Considering the locations of SKI and 
CHH 12/13, however, HL 8 could be the counterpart of the CHH 
12-13 outflow. 

HL 9.— At -8" south of IRAS 4BI, the two peaks of HL 
9 line up in the north-south direction, which agrees with the 
direction of the IRAS 4BI bipolar outflow seen in the HCN line 
(Choi 2001b). 

ASR 57. — This is one of the brightest H2 emission features in 
the region imaged. The location of ASR 57 at the intersection of 
the IRAS 4A outflow and the HH 7-1 1-HL 7 outflow raised the 
possibility that the two flows might be interacting (Hodapp & 
Ladd 1995). Since the dense core obstructing the IRAS 4A out- 
flow is a part of the large cloud associated with SVS 13 and the 
HH 7-1 1 complex (Paper I; Choi 2005b), the two outflows may 
be close to each other along the line of sight. Further studies 
are needed to test whether or not they are physically interacting 
indeed. Proper motion studies of the H2 peaks within ASR 57 
may be useful. 

VIA 13/30/32. — These centimeter continuum sources were 
detected in the UH image as point-like sources. The agreement 
between the VLA position and the near-IR position is good 
within 0'.'2. 



4.5. Other H2 Emission Features 

The H2 images (Figs. 1 and 2) show several emission fea- 
tures that were previously undetected or unlabeled. Coordinates 
of their peak positions are listed in Table 1 . These H2 features 
are labeled with a prefix CHH and briefly described here. Some 
of the previously known H2 features are also discussed below. 

CHH 1. — This is a short jet-like H2 feature. A probable driv- 
ing source is VLA 14. 

CHH 2/4.— These H2 features and HH 347B probably be- 
long to a single outflow. 

CHH 3. — This is a short jet-like H2 feature. A probable driv- 
ing source is ASR 74, the star just north of CHH 3. 

CHH 5. — This is a linear chain of H2 emission features. Ho- 
dapp & Ladd (1995) suggested that CHH 5 and HL 4 are re- 
lated with the SVS 13 complex. Another possibility is that the 
CHH 5-HL 4 flow may be driven by VLA 19 (SK 14). Yet an- 
other possibility is that CHH 5b, which is elongated roughly in 
the northeast-southwest direction, may be a short outflow driven 
by VLA 19, while CHH 5a/c (and HL 4) may belong to a sepa- 
rate outflow. 

CHH 6. — This is an unusually straight jet-like H2 feature. It 
appears to be related with ASR 20/21 and probably with HH 
344. 
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